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Abstract. In this paper we consider the effect of a self-consistently computed magne- 
tosdisc field structure on the magnetosphere-ionosphere coupling current system at Jupiter. 
We find that the azimuthal current intensity, and thus the stretching of the magnetic 
field lines, is dependent on the magnetosphere-ionosphere coupling current system pa- 
rameters, i.e. the ionospheric Pedersen conductivity and iogenic plasma mass outflow rate. 
Overall, however, the equatorial magnetic field profiles obtained are similar in the in- 
ner region to those used previously, such that the currents are of the same order as pre- 
vious solutions obtained using a fixed empirical equatorial field strength model, although 
the outer fringing field of the current disc acts to reverse the field-aligned current in the 
outer region. We also find that, while the azimuthal current in the inner region is dom- 
inated by hot plasma pressure, as is generally held to be the case at Jupiter, the use of 
a realistic plasma angular velocity profile actually results in the centrifugal current be- 
coming dominant in the outer magnetosphere. In addition, despite the dependence of the 
intensity of the azimuthal current on the magnetosphere-ionosphere coupling current sys- 
tem parameters, the location of the peak field-aligned current in the equatorial plane also 
varies, such that the ionospheric location remains roughly constant. It is thus found that 
significant changes to the mass density of the iogenic plasma disc are required to explain 
the variation in the main oval location observed using the Hubble Space Telescope. 



1. Introduction 

The dynamics of Jupiter's middle magnetosphere 
are dominated by planetary rotation coupled with the 
centrifugally-driven outflow of plasma from the volcanic 
moon lo, which orbits at ~5.9 Rj (where Rj is the equa- 
torial radius of Jupiter equal to 71,373 km) and liberates 
~1000 kg s~^ of sulphur and oxygen into a torus surround- 
ing the satellite's orbit [e.g. Siscoe and Summers, 1981; Va- 
syliunas, 1983; Khurana and Kivelson, 1993; Delamere and 
Bagenal, 2003; Dols et ai, 2008]. Subrotation of outflow- 
ing equatorial plasma leads to the bend-back of magnetic 
fleld lines out of the meridian planes and the formation 
of a large-scale magnetosphere-ionosphere (M-I) coupling 
current system, illustrated schematically in Figure 1 [e.g. 
mil, 1979, 2001; Pontius, 1997; Cowley and Bunce, 2001; 
Nichols and Cowley, 2003, 2004, 2005]. The current sys- 
tem, which consists of an equatorward-directed Pedersen 
current in the ionosphere and a radial current in the equa- 
torial plane joined in the inner region by an upward-directed 
fleld-aligned (Birkeland) current and closed in the outer re- 
gion by a downward-directed current, communicates drag 
from the atmospheric neutrals to the equatorial plasma. The 
upward fleld-aligned component of this current system, as- 
sociated with downward-precipitating electrons, is thought 
to be the cause of Jupiter's main auroral oval, which is the 
most signiflcant of Jupiter's various ultraviolet (UV) auro- 
ral forms [Grodent et ai, 2003a; Clarke et al, 2004; Nichols 
et ai, 2009]. 

This jovian M-I coupling current system was studied orig- 
inally by Hill [1979], who calculated the plasma angular ve- 
locity proflle employing a dipole planetary magnetic fleld. 
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and the theory was later generalised to include a realistic 
current sheet magnetic field (magnetodisc) model by Pon- 
tius [1997]. The link with the main oval auroral emission 
was realised later by Hill [2001] , who again used a theoretical 
angular velocity profile calculated using a dipole field, and 
Cowley and Bunce [2001], who used empirical plasma angu- 
lar velocity and current sheet magnetic field profiles. Cowley 
et al. [2002] computed using Hill-Pontius theory the plasma 
angular velocity and current profiles using both dipole and 
current sheet field models, and showed that the stretching of 
the equatorial middle magnetosphere field lines associated 
with the current sheet dramatically alters the magnitude 
and location of the auroral field-aligned currents. The ef- 
fects of two poorly-constrained system parameters, the effec- 
tive ionospheric Pedersen conductance Ep, and the plasma 
mass outfiow rate M were then studied in detail by Nichols 
and Cowley [2003], and Nichols and Cowley [2004] went on 
to examine the effect of self-consistent modulation of the 
ionospheric Pedersen conductance due to auroral electron 
precipitation. Nichols and Cowley [2005] and Ray et al. 
[2010] have since studied the effect of field-aligned voltages, 
and Cowley et al. [2007] examined the effects of solar wind- 
induced expansions and contractions of the planet's magne- 
tosphere. In addition, the modulation of the current system 
by diurnal variation of the ionospheric Pedersen conduc- 
tance caused by solar illumination has been studied by Tao 
et al. [2010]. Most recently, Nichols [2011] applied approxi- 
mations derived by Nichols and Cowley [2003] for the jovian 
M-I coupling current system to the cases of rapidly rotat- 
ing, strongly illuminated Jupiter-like exoplanets, the radio 
emissions from which may offer a novel method for detecting 
such objects. 

A key limitation of the preceding studies, however, is that 
they have all employed fixed magnetic field models as the 
basis for the computations, be they dipole or current sheet 
in form. However, as shown in Figure 2 (reproduced from 



1 



X - 2 



J. D. NICHOLS: JOVIAN M-I COUPLING AND MAGNETODISC 




Figure 1. Sketch of a meridian cross section through 
Jupiter's inner and middle magnetosphere, showing the 
principal physical features involved. The arrowed solid 
lines indicate magnetic field lines, the arrowed dashed 
lines the magnetosphere-ionosphere coupling current sys- 
tem, and the dotted region the rotating disc of outflowing 
plasma. After Cowley and Bunce [2001]. 

Grodent et al. [2008]), the main oval has been observed to 
shift in latitude by up to ~3° when comparing images ob- 
tained a number of years apart. This shift in the latitude of 
the main oval was accompanied by a similar shift of ~2° in 
the latitude of the footprint of Ganymede, such that Grodent 
et al. [2008] attributed the shift in the main oval to a change 
in the intensity of the azimuthal current, which affects the 
mapping between ionosphere and equator, rather than sim- 
ply a shift across L-shells of the field-aligned current, as is 
predicted would occur for different values of Ep and M [e.g. 
Nichols and Gowley, 2003]. Caudal [1986] showed that the 
stretching of Jupiter's middle magnetosphere field lines is 
caused in part by the centrifugal force of iogenic plasma, a 
quantity specifically associated with the iogenic plasma mass 
density and angular velocity profile. He constructed a model 
for Jupiter's magnetodisc by modifying a terrestrial storm- 
time ring current model [Lackner, 1970], and employing in- 
puts based on Voyager plasma temperature and density ob- 
servations, along with a fixed angular velocity profile given 
by the theory of Hill [1979], computed using a dipole field 
model. Caudal [1986] did note the inconsistency in employ- 
ing an angular velocity profile calculated using a dipole field, 
but while Pontius [1997] and Cowley et al. [2002] showed 
that the equatorial plasma angular velocity profile is rela- 
tively insensitive to the field model used, the latter authors 
showed that the resulting auroral currents are very sensi- 
tive to the model employed. In addition, in calculating the 
iogenic plasma angular velocity profile Caudal [1986] used 
a corotation breakdown scale distance pH (termed the 'Hill 
distance') of 20 Rj in conformity with the value deduced 
by Hill [1980], and in his model the plasma angular velocity 
thus falls to ~17% of rigid corotation by 60 Rj. However, ob- 
servational studies such as Kane et al. [1995] have reported 
that the plasma angular velocities remain at ~50% of rigid 
corotation out to ~60 Rj, and Hill [2001] later revised his es- 
timate of pH to 30 Rj, such that the middle magnetosphere 
plasma angular velocities employed by Caudal [1986], and 
thus the centrifugal force imparted by the iogenic plasma, 
are somewhat lower than realistically expected. The purpose 
of the present paper is thus as follows. First, we incorpo- 
rate the calculation of the plasma angular velocity profile 
using Hill-Pontius theory into the model of Caudal [1986], 
such that the resulting magnetosphere-ionosphere currents 
are computed using values of the equatorial magnetic field 
self-consistent with the plasma angular velocity profile. Sec- 
ond, in doing so we will update the model results of Cau- 
dal [1986] using more realistic plasma parameters, including 
values obtained from Galileo data. We then examine the 
effect of the ionospheric Pedersen conductance and iogenic 



plasma mass outflow rate in order to compare results with 
the previous modeling work of Nichols and Cowley [2003], 
and examine whether variations of these parameters may be 
responsible for the changing auroral locations observed by 
Grodent et al. [2008]. 

2. Theoretical background 

2.1. M-I coupling current system equations 

We begin by discussing the equations governing the jo- 
vian M-I coupling current system. The system has been dis- 
cussed in depth previously [e.g. Hill, 1979; Pontius, 1997; 
Cowley et al, 2002; Nichols and Cowley, 2003, 2004, 2005], 
such that only a brief outline is given here. We flrst assume 
axisymmetry, such that Jupiter's poloidal field can be de- 
scribed with the use of a flux function F{p,z), related to 
the Euler potential a used in Section 2.2, and from which 
the magnetic field can be computed via 




Figure 2. Top panel: Superposition of the polar pro- 
jection of two images of Jupiter's northern aurora ob- 
tained with HST more than four years apart. The red 
image was obtained with the Space Telescope Imaging 
Spectrograph (STIS) in December 2000, and the blue im- 
age was obtained with the Advanced Camera for Surveys 
(ACS) camera in AprU 2005. The 90° and 180° System 
HI meridians have been highlighted on a 10 spaced grid. 
Green arrows point to the footprints of Ganymede and 
lo, and the main emission has also been marked in green. 
Bottom panels: Individual polar projections using the 
same longitude system as in top panel. From Grodent 
et al. [2008]. 
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B=Q)vFx^, (1) 

where p is the perpendicular distance from the magnetic 
axis, 2 is the distance along this axis from the magnetic 
equator, and is the azimuthal angle. Mapping between 
the equator and ionosphere is then easily achieved by writ- 
ing Fe = Fi, where subscript 'e' refers to the equator and 'i' 
refers to the ionosphere. Although the magnetic field model 
discussed in Section 2.2 does compute the field at all lat- 
itudes and altitudes down to 1 Rj, such that in principle 
it provides ionospheric field values, Caudal [1986] pointed 
out that the simple treatment of the magnetopause currents 
renders the model invalid at latitudes above ~ 50°, and in 
addition the cartesian grid used here is too coarse to be 
useful in this region. However, the ionospheric field is over- 
whelmingly dominated by the planetary dipole, such that 
in common with previous works we assume the field in the 
ionosphere is purely dipolar. The ionospheric flux function 
is then given by 

Fi=Bjp1 = BjR'j8in^ei , (2) 

where Bj is the dipole equatorial magnetic field strength 
(equal to 426, 400 nT in conformity with the VIP 4 internal 
field model of Connerney et al. [1998]), Rj is Jupiter's radius 
(equal to 71 373 km), pi is the perpendicular distance from 
the magnetic axis, and 0i represents magnetic co-latitude. 

The application of Newton's second law to an axisymmet- 
ric radially-outward steady flow of plasma from the torus 
source yields the 'Hill-Pontius' differential equation for the 
plasma angular velocity ui, given by 



Recalling that Fi = Fe, we have from Eqs. 4, 5 and 2 

such that the total radial current integrated in azimuth Ip 

is 

Ip = 27TpJp = 87rSpnjF, (^1 - , (7) 

which is equal to twice the azimuth-integrated equatorward- 
directed Pedersen current Ip flowing in each hemisphere. 
The field-aligned current density at the top of the iono- 
sphere is then computed from the divergence of either 
total field-perpendicular current, such that, in terms of the 
radial current 

• — dip , , 

2.2. Magnetodisc field model 

Whereas previous studies have specified the equatorial 
magnetic field profile as a fixed input to the equations dis- 
cussed in Section 2.1, here we self-consistently employ the 
jovian magnetodisc model of Caudal [1986], which we note 
has also recently been adapted to the saturnian magnetodisc 
[Achilleos et al., 2010]. Briefly, Caudal's [1986] model rep- 
resents the spin-aligned magnetic field as the gradient of the 
Euler potentials, given generally by a and 13, but which in 
the axisymmetric case are reduced to one function a{r,6), 
where (r, 6) are polar coordinates. The Euler potential a is 
related to the flux function F via 
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where pe represents equatorial radial distance, f2j is the 
planet's angular velocity equal to 1.76 x 10^* rad s~^, and 
\Bze\ is the magnitude of the north-south magnetic field 
threading the equatorial plane. Note that the effective Ped- 
ersen conductance Ep (here defined for one hemisphere) is 
reduced from the true value Ep by Sp = (1 — fc)Ep, where 
the parameter k represents the reduction of the angular ve- 
locity of the neutral atmosphere (O}) from rigid corotation 
(fij) via 'slippage' [Huang and Hill, 1989; Millward et al., 
2005], such that (fij - fi}) = k{Qj - oj). The value of k is 
somewhat uncertain, so in common with previous studies we 
take k = 0.5, although we note that in reality this approach 
may be an oversimplification [Smith and Aylward, 2009; Tao 
et al., 2009]. The quantities Fe and \Bze\ are obtained from 
the magnetodisc model discussed in Section 2.2, such that 
Eq. 3 is solved numerically to obtain the equatorial plasma 
angular velocity profile. 

We now discuss the equations which describe the result- 
ing magnetosphere- ionosphere coupling currents. First, the 
equatorward-directed height-integrated Pedersen current ip 
is given by 



ip = 2T.pBjQ.jpi 1 



Qj 



(4) 



where we have taken the ionospheric field to be vertical and 
equal to 2Bj in strength (an approximation valid to within 
~5% in our region of interest [Nichols and Cowley, 2003]). 
Current continuity and the assumption of north-south sym- 
metry then yields for the equatorial radial current integrated 
across the width of the current sheet io 



peip = 2piip 



(5) 



(9) 



such that the magnetic field components axe then given e.g. 
in cylindrical coordinates by 



Rj da 
p dz 
Rj da 
p dp 

B^ = . 



Bo = 



B, = 



(10a) 

(10b) 

(10c) 



Note that this assumes that the poloidal magnetodisc struc- 
ture is unaffected by the small azimuthal field generated by 
the equatorial radial component of the M-I coupling current 
system. Caudal [1986] considered the momentum equation 
for a rotating plasma, i.e. 

j X B = VP-dui^pp , (11) 

where j is the current density, P is the plasma pressure and 
d is the plasma mass density. Prom this he derived the dif- 
ferential equation 



(1 



= -9{r,x,a) , 



(12) 



where x = cos (9 and the function g, which is derived from 
the plasma pressure and angular velocity distributions, is 
related to the azimuthal current density via 

g = po3vP ■ (13) 

Function g comprises two summed components, represent- 
ing contributions from the hot ~30 keV [Krimigis et al., 
1981] and cold ~100 eV [McNutt et al, 1981; Frank et al, 
2002] plasma populations in Jupiter's magnetosphere. For 
the hot plasma population (subscript 'h') the pressure grar 
dient dominates the centrifugal force, such that the latter is 
neglected and g is given by 
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9h = Ho 



R-j) 



(14) 



while the source function for the cold population (subscript 
'c') includes the centrifugal force, such that 



■ Mo 



P 

Rj 



cxp 



2/2 



dPc o ^ PcoRj 

da 



P\B, 



(15) 



where po, Pco and Bzeo arc the values at the equatorial 
crossing point of the field line. Quantity I in Eq. 15 repre- 
sents the centrifugal equatorial confinement scale height of 
the cold plasma, given for a singly-ionised, monoionic pop- 
ulation with temperature Tc and ion mass m by 



^ /' 2kTc 



(16) 



where k is Boltzmann's constant equal to 1.38x10"^"^ J K~^. 
The plasma pressure P is given in general by the ideal gas 
law for a singly- ionised plasma 



employed Hill's [1979] solution to Eq. 3 for the dipole field, 
given by 



n7; = ^^^P 
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(23) 



where erf (z) = (2/-^/7r) e * dt is the error function and 
the Hill distance ph is given by 



m 



f 2-KT.''pB/R/ 
V M 



(24) 



which, as discussed in Section 1, Caudal [1986] took 
to be equal to 20, corresponding to a value for the 
quotient (Ep/M) = 2.75 x 10"^ mho s kg"^ (where 
1 mho = 1 siemcn). It is also worth noting here that 
Hill's [2001] revised estimate of [pn/Rj) = 30 corresponds 
to (Ep/M) = 1.4 X 10"" mho s kg-^ 



2NkT 
V 



(17) 



where A'^ is the number of ions per Wb and V{a) is the 
volume of the unit flux tube, given for the hot plasma by 



Vh 



-I- 



and for the centrifugally confined cold plasma by 



\4 = 



cxp 



2/2 



ds 



(18) 



(19) 



where the exponential term in equation 19 represents the 
centrifugal confinement pressure of the cold plasma. All the 
physical properties of the plasma are thus represented in P 
and I, and Caudal [1986] used Voyager observations [Bage- 
nal and Sullivan, 1981; Connemey et al., 1981; Krimigis 
et al., 1981; McNutt et al., 1981; Siscoe and Summers, 1981] 
to provide suitable values. Specifically, he took for the hot 
plasma 



{_ 3.0 X 10'^ 
Vh{a) 
OC Po 



if Po > 7.5 Rj 
if Po < 7.5 Rj , 



(20) 



a form which we also employ here. For the cold plasma he 
employed equations 16, 17, and 19, with profiles for Nc{po) 
and fcTc(po) derived from Voyager data, i.e. 



iVc(po) = < 



if Po < 5 Rj 

10.7 X 10^2 if 5.7 Rj < Po < 7 Rj (21) 

2.9 X 10^^ if Po > 8 Rj , 



with continuity achieved through linear interpolation be- 
tween these domains, and 



fcrc(po) 



1 eV if Po = 5 Rj 

35 eV if 6 Rj < Po < 7 Rj (22) 

10{po/Rj) eV if Po > 9 Rj , 



with continuity achieved here through linear interpolation 
of log(fcTc). For the plasma angular velocity. Caudal [1986] 



The iterative analytic solution to equation 12, stated 
by Caudal [1986] and derived explicitly by Achilleos et al. 
[2010], is 



On = CKo + (1 



^ 2n-|-3 



71 = 

-^(^|^\"+V(«)d«)+r"+^£ 



9n{u)u 



(25) 



where Pi'^'''\z) are the Jacobi polynomials [see e.g. 
Abramowitz and Stegun, 1965], and 



1 

9n{u) = -f^ J ^g{r,x)P2'^\x)dx 



where 



{l-x^){Pl,''^\x)fdx 



(26) 



(27) 



The solution is initiated with the Euler potential for a dipole 
field ao given by 



oio = BjRj ' 



(28) 



and proceeds by iteration according to the scheme illustrated 
by Figure 4 in Caudal [1986]. As he noted, in order to 
achieve convergence, after a few iterations the new values of 
a are obtained using a weighted average of a„ and a„_i. 
In addition, at each iteration an Euler potential repre- 
senting the field induced by the equatorial magnetopause 
current as seen inside the magnetosphere is added to the so- 
lution given by equation 25. It is modeled as an irrotational 
field of strength Bs , such that 



and Bs is given by 



2Rjamp 



: 0.6 



R'^ 



(29) 



(30) 



where Rmp is the distance to the equatorial magnetopause, 
taken in this study to be the representative value of 85 Rj, 
and Ump is the value of a at the equatorial magnetopause. 
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Figure 3. Plot showing tiie number of cold ions per We- 
ber versus radial distance. The solid line shows the values 
employed in this study, the dotted line shows the values 
used by Caudal [1986], and the dashed line shows the 
values estimated using the results of Frank et al. [2002] 
along with the flux tube volume calculated using Cau- 
dal' s [1986] model. 



2.3. Application of the magnetodisc model to M-I 
coupling 

In applying Caudal's [1986] magnetodisc model to the jo- 
vian M-I coupling current system, the major development 
of the model is the treatment of the cold plasma angular 
velocity, which we describe below. We first discuss, how- 
ever, a secondary modification of the model, concerning the 
input cold plasma parameter values. As discussed in Sec- 
tion 2.2, Caudal [1986] employed cold plasma ion number 
density values (number per Wcbor) based on estimates over 
the radial range 5-9 R,i calculated by Bagenal and SulU- 
van [1981], who used Voyager plasma data and assumed a 
dipole field and an exponential distribution of plasma along 
the field lines. Between 8 Rj and the magnetopause Caudal 
[1986] used a constant value based on the outer values of the 
Bagenal and Sullivan [1981] results. More recently, Frank 
et al. [2002] have reported thermal plasma density observa- 
tions obtained by the Galileo spacecraft over a much greater 
radial distance, out to 100 Rj, and provided the following 
power laws for the thermal plasma number densities 



nc(po) 



3.2 X 10^ po" 
9.8 po"^-^* cm-3 



cm-3 if po < 20 Rj 
if Po > 50 Rj 



(31) 



Estimates of the number of ions per Weber can then be ob- 
tained by multiplying equation 31 by the weighted flux tube 
volume (e.g. as given by equation 19), for which the values as 
calculated using Caudal's [1986] original model may be used 
as reasonable estimates. The number of ions per Weber thus 
calculated are shown in Figure 3. The dotted line (which 
is essentially overlaid by the solid line within 8 Rj) shows 
the values employed by Caudal [1986], which we recall are 
unconstrained by data beyond 9 Rj, while the dashed line 
shows the estimates using the two power laws given by Frank 
et al. [2002] as discussed above, where wc switch from one to 
the other at their intersection at ^^22 Rj. It is apparent that, 
although the estimates using the Frank et al. [2002] profile 
is in reasonable agrccinciit at ~8 Rj with the Bagenal and 
Sullivan [1981] results, the constant value taken by Caudal 
[1986] in the region beyond significantly overestimates the 
Frank et al. [2002] results over most of this region. We thus 
employ a revised estimate of the number of ions per Weber 
beyond 8 Rj given by the mean of the log of the Frank et al. 
[2002] values, i.e. 10<'°s^=> = 8.1 x 10^^ ions per Weber. 



A constant value has been maintained both for simplicity 
and since, while the flux tube volumes obtained by Caudal 
[1986] are the most reasonable to use, the values have not 
been verifled experimentally, such that here we simply use 
them to obtain an appropriate spot value for the number 
of ions per Weber. The values of Nc used in this study 
are therefore shown by the solid line in Figure 3. A second 
minor modification of the thermal plasma parameters con- 
cerns the plasma temperature. As shown in equation 22, 
Caudal [1986] used fcTe(po) = 10{pa/Rj) eV beyond 9 Rj. 
However, Frank et al. [2002] report thermal plasma temper- 
atures of ~500 eV at ~25 Rj, such that here wc instead 
employ kTc{pa) = 20(po/i?j) eV beyond 9 Rj. Overall, 
therefore, on the basis of the Galileo results presented by 
Frank et al. [2002], the cold plasma in the model used in 
this study is slightly warmer and less dense than that used 
in Caudal's [1986] model. 

We now discuss the treatment of the plasma angular ve- 
locity. As discussed in Section 2.2, Caudal [1986] used the 
fixed angular velocity profile for a dipole field given by equa- 
tion 23 at every iteration of the model. Here, wc also initiate 
the model with a dipole magnetic field across the entire do- 
main and then iteration proceeds according to his Figure 4, 
except that we instead solve radially outward numerically 
at every iteration equation 3 with F^ and \Bze\ computed 
from the equatorial values of a„-i using equations 9 and 
10b, respectively, to obtain the angular velocity profile con- 
sistent with that iteration of the magnetic field, rather than 
simply always using the dipole solution of the Hill-Pontius 
equation. The resulting angular velocity profile is then em- 
ployed in equation 16, and thence the function gc given by 
equation 15. Thus, the current sheet magnetic field and an- 
gular velocity profiles are always consistent with each other, 
and evolve together toward a converged self- consistent solu- 
tion. The resulting M-I currents are calculated upon con- 
vergence, defined as being when the max:imum relative dif- 
ference in the values of a between one iteration and the next 
is less than 0.5%. In some cases, the extra degree of free- 
dom introduced by the modification of the plasma angular 
velocity profile between successive iterations results in the 
model reaching a 'quasi-steady' state, rather than true con- 
vergence, in which the model perpetually fluctuates about 
a set of values. In these cases the model run is stopped and 
the mean of the results of the last 5 iterations is used, a 
number which yields results representative of the set of pro- 
files obtained when the model has reached a quasi-steady 
state. 

In the results that follow we have employed ranges of val- 
ues of the ionospheric Pcdersen conductance and the iogenic 
plasma mass outflow rate, two parameters whose exact val- 
ues arc unknown and the effects of which on the M-I coupling 
current system (assuming a fixed magnetic field model) have 
been studied previously [Nichols and Cowley, 2003]. In this 
study, the Pedersen conductance is assumed for simplicity to 
be constant, since this allows easy comparison with the anar 
lytic results of Nichols and Cowley [2003], although we note 
that in reality feedback resulting from the precipitating elec- 
tron flux will enhance the conductivity in the auroral region, 
modifying the plasma flow and current profiles as shown by 
Nichols and Cowley [2004]. The iogenic mass outfiow rate 
can be treated in two ways. First, if the cold ion number 
density is assumed constant, the mass outflow rate is then 
related solely to the rate of outward transport of iogenic 
plasma, such that higher mass outflow rates equate to faster 
outward transport. On the other hand, if the outward trans- 
port rate is instead assumed constant, the cold ion number 
density is proportional to the mass outflow rate. Taking the 
canonical value of M = 1000 kg s~^ as a reference, the cold 
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Figure 4. Plot showing tlie magnetic field and cur- 
rent sheet structures as computed using 3 values of the 
quotient (EJ, / M), i.e. (a) 10-^ (b) 10"'', and (c) 
5 X 10~* mhoskg"^. The black lines are contours of 
a, thus indicating magnetic field lines, and the colours 
indicate the azimuthal current density in uA m""^. 



ion number density values described above are then modified 
by 



M 



(32) 



1000 kg s- 

In Section 3 we thus show results for both these scenarios, 
and we note that the reality will probably lie somewhere 
between these two cases. 

3. Results 

3.1. Results with N^. independent of M 

We now present the results obtained using the model de- 
scribed in Section 2. We first show in Figure 4 the structures 
of the magnetic field (black contours) and total azimuthal 
current density (colours) as computed using 3 values for the 
quotient (SJ,/M) = 10"^, IQ-*, and 5 x 10"" mhoskg"^ 
These roughly bracket both the value assumed by Caudal 
[1986] and the revised value of Hill [2001] as discussed in 
Section 2. Note that here we keep the cold plasma den- 
sity independent of the plasma mass outflow rate. Hill 
[1979] showed that higher values of (Ep/M) result in higher 
plasma angular velocity values, and it is apparent that 
higher values of this quotient result in a more stretched 
magnetic field structure with a thinner, more intense cur- 
rent sheet, particularly in the region outward of ~20 Rj. 
Specifically, the half-width of the current sheet in the mid- 
dle magnetosphere is typically ~8-10, ~6-8, and ~3-5 Rj 



for (Ep/M) = 10" 



10" 



and 5 X 10 mho s kg , re- 



spectively. These are all somewhat larger than the value 
of 2.5 Rj employed in the empirical 'Voyager-l/Pioneer-10' 
('CAN') current sheet field model of Connerney et al. [1981], 
with the (Ep/M) = 5 x 10~" mho s kg~^ result being most 
consistent with the latter. The increased azimuthal current 
for higher values of (Ep/M) is required to balance the ele- 
vated centrifugal force imparted by the faster-rotating equa- 
torial plasma for higher values of (Ep/Af). 

This can be further appreciated from Figure 5, in which 
we show various parameters associated with the magne- 
todisc model and M-I coupling current system for each of 



the above values of (Ep/M), and where for the purposes of 
the M-I coupling current calculations we take the canoni- 
cal value of M = 1000 kg 8"^ Specifically, we show, from 
top to bottom, the magnitude of the north-south magnetic 
field threading the equatorial plane |_Bze| in uT, the iono- 
spheric co-latitude to which the magnetic field maps 9i in 
degrees, the equatorial plasma angular velocity normalised 
to the planet's rotation rate (u/Qj), the ratio of the equa- 
torial azimuthal current density associated with the cold 
plasma centrifugal force to that of the hot plasma pressure 
(ic/)o cent/jVo h), the cold plasma pressure Pco in Pa, the 
azimuthally-integrated equatorial radial current Ip in MA, 
and finally the field-aligned current density jyj at the top 
of the ionosphere in /.lA m~^, all versus equatorial radial 
distance in Rj. Note that the solid coloured lines indicate 
results from model runs which converged, while long-dashed 
lines indicate results from model runs which have reached 
a 'quasi-steady' state as discussed in Section 2.3. Starting 
with the equatorial magnetic field strength |-Bze| shown in 
Figure 5a, it is first evident that all three model results are 
similar out to ~15 Rj, beyond which they diverge. Also 
shown in Figure 5a for comparison are the magnetic field 
strengths given by the pure planetary dipole (dashed black 
line), given by 



Hze dip {Pe 



Rj_ 

pe 



(33) 



and the 'CAN-KK' current sheet magnetic field model of 
Nichols and Cowley [2004] (dot-dashed line), given by 



B 



ze CAN-KK 



B'{^ 



exp 



5/2' 



pe 



(34) 



where B'^ = 3.335 x 10^ nT, p* = 14.501 Rp, Bo = 
5.4 X 10* nT, and m = 2.71. This form closely approx- 
imates the field model used by Cowley and Bunce [2001] 
and Cowley et al. [2002, 2003], who employed the CAN field 
model of Connerney et al. [1981] in the inner region and 
the Voyager-1 ('KK') outbound pass model of Khurana and 
Kivelson [1993] in the outer region. In the model results 
obtained in this study, the |-Bze| values are less than those 
for the dipole in the inner region owing to the radial dis- 
tention of the field by the current sheet. All three results 
are reasonably consistent with the CAN-KK model to dis- 
tances of ~20 Rj, but (Ep/M) = 5 X 10"" mhoskg"^ 
again gives the best agreement, roughly tracking the CAN- 
KK model values out to ~40 Rj. We note that the slight 
jitter in the latter \Bze\ profile between ~40 and ~60 Rj 
is representative of the spontaneous instability which pro- 
hibits models runs with higher values of (Ep/Af) from truly 
converging, and which, as mentioned by Caudal [1986], can 
lead to the formation of neutral points for more stretched 
magnetodiscs. In the outer region, the \Bze\ values become 
greater than those for the dipole, with the transitions occur- 
ring at ~38, ~42, and ~52 Rj for {T,*p/M) = 10-^ 10"", 
and 5 x 10~" mho s kg"'^, respectively. This transition, also 
originally noted by Caudal [1986], is due to the outer fring- 
ing fields of the current sheet, and the outer- most values of 
l-Bzej of ~10-20 nT are consistent with the values of ~16 nT 
observed just inside the magnetopause [Acuna et al., 1983]. 
Note that in contrast, in the radial range of Figure 5 the 
'CAN-KK' values are always less than those of the dipole, 
indicating that the current sheet in the pre-dawn region of 
the Voyager-1 outbound pass was evidently extended due to 
the distant (~160 Rj) magnetopause in this region [Acuna 
et al, 1983]. 
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Figure 5. Plot showing magnetodisc and current system 
parameters computed for (Ep/M) — 10~^ (blue), 10~* 
(green), and 5 x 10~* (red) mho s kg~^, plotted versus 
equatorial radial distance. Parameters shown are (a) the 
magnitude of the north-south magnetic field threading 
the equatorial plane |i3ze| in nT, (b) the ionospheric co- 
latitude to which the magnetic field maps 6i in degrees, 
(c) the equatorial plasma angular velocity normalised to 
the planet's rotation rate (uu/Qj), (d) the ratio of the 
equatorial azimuthal current density associated with the 
cold plasma centrifugal force to that of the hot plasma 
pressure (i^o cent/j</>o h), (e) the cold plasma pressure 
Pco in Pa, (f) the azimuthally-integrated equatorial ra- 
dial current Ip in MA, and (g) the field-aligned current 
density at the top of the ionosphere in fj,A m~^. In 
Figures 5a and 5b the dashed black lines show the plan- 
etary dipole values, and the dot-dashed black lines show 
the values using the model of Nichols and Cowley [2004]. 
In Figure 5c the horizontal dotted line indicates rigid 
corotation and the dashed black line indicates the profile 
taken by Caudal [1986]. In Figure 5d the horizontal dot- 
ted line indicates where the azimuthal current densities 
associated with the cold plasma centrifugal force and the 
hot plasma pressure are equal. In Figure 5e the dashed 
and dot-dashed black lines show the power laws given by 
Frank et al. [2002] . Note that both the positive and nega- 
tive values in Figure 5g are plotted on logarithmic scales, 
such that the horizontal line is at ±0.0001, and the re- 
sulting apparent discontinuities at the transition points 
are simply artefacts of the plotting scale. In Figures 5f 
and 5g the canonical value of M = 1000 kg s~^ is used. 
In all panels the solid coloured lines indicate results from 
model runs which converged, while long-dashed lines indi- 
cate results from model runs which have reached a 'quasi- 
steady' state as discussed in Section 2.3. In Figures 5c, 
5f and 5g the coloured dot-dashed lines indicate results 
obtained using the fixed magnetic field model of Ntchols 
and Cowley [2004] for comparison. 



Figure 5b shows the ionospheric co-latitude to which the 
magnetic field maps, calculated using equations 2 and 9. 
Also shown are the values for the planetary dipole (dashed 
black line), for which 



i^eCAN-KK(pe) = 



2.5p- 



5/2' 



+ 



Bo 



R., 



(m - 2) \ pe 



(36) 



Fedip{pe) 



BrRi 



(35) 



where Fao ~ 2.841 x 10'' nT Rj^ is the value at infinity, and 
r(a, z) = f"~^e~' dt is the incomplete gamma function. 
It is apparent that for each value of (Ep/A/) used the field 
line mapping is more consistent with that of the CAN-KK 
field model than the dipole, although the elevated values of 
\Bzs\ in the outer region relative to the CAN-KK values re- 
sults in a broadening of the ionospheric latitudinal band to 
which the outer magnetosphere maps. It is evident, however, 
that for increased values of (Ep/M) the middle magneto- 
sphere field lines map to a modestly more equatorward and 
thinner latitudinal band in the ionosphere. For example, 
field lines threading the equatorial plane between 20-60 Rj 
map to between ~12.6-16.6°, ~14.0-16.9°, ~15. 6-17.1° for 
(Ep/M) = 10"^ 10"*, and 5 x 10"" mho s kg-\ respec- 
tively, with the latter result being most consistent with the 
CAN-KK model. In addition, it is worth noting that in this 
model the ionospheric co- latitudes of the last closed field line 
are ~7.9°, ~9.0°, and ~11.0° for (Ep/M) = 10"^ 10"*, 
and 5 x 10"* mho s kg"^, respectively. The latter value is 
in excellent agreement with the value of ~11° recently deter- 
mined by Vogt et al. [2011], and is also consistent with the 
value of 10.25° used by Cowley et al. [2005] in their global 
model of Jupiter's polar ionospheric flows. 

The angular velocity of the equatorial plasma as com- 
puted by the present model is shown in Figure 5c, along 
with the angular velocity profiles calculated using the fixed 
CAN-KK field model and the profile assumed by Caudal 
[1986] for comparison. In all cases the profiles calculated 
here are similar to those obtained using the CAN-KK field 
model in the inner region, and deviate toward higher values 
in the outer region due to the increased j x B force owing to 
the elevated values of [Bzej relative to the CAN-KK model. 
As mentioned in Section 2, the angular velocity profile used 
by Caudal [1986] is roughly consistent with the lowest value 
of (Ep/M) used here, and reduces to ~0.17 by 60 Rj. On 
the other hand, observational studies such as Kane et al. 
[1995] have reported that the plasma angular velocities re- 
main at ~0.5 out to ~60 Rj. The present angular velocity 
profile which best fits this behaviour is that produced using 
(Ep/Af) = 10"* mho s kg"'^, which is also generally consis- 
tent with the values at the plasma sheet crossings obtained 
by McNutt et al. [1981] (i.e. those values at the local max- 
ima in the Voyager 1 data shown in their Figure 21). Lower 
and higher values of (Ep/Af) then produce angular velocity 
profiles which are overall respectively somewhat lower and 
higher than observations suggest. 

The effect of the plasma angular velocity on the azimuthal 
current is shown in Figure 5d, in which we plot the ra- 
tio of the equatorial azimuthal current density associated 
with the centrifugal force to that of the hot plasma pressure 
( ^0 o cent /i</) o h) , giving an indication as to which of these two 
components of the azimuthal current is dominant. Caudal 
[1986] concluded that the latitude-integrated current associ- 
ated with the hot plasma pressure dominates both the cold 
plasma pressure current and the centrifugal force current 
over the whole of the magnetosphere. This result was sup- 
ported by Achilleos et al. [2010], although these authors also 
showed that in the original Caudal [1986] model the effect of 
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the centrifugal force strongly peaks near ~27 Rj, such that 
equatorial current densities associated with the hot plasma 
pressure and centrifugal force become comparable between 
~20-30 Rj, a concern which was originally raised by Mauk 
and Krimigis [1987] on the basis that it apparently contra- 
dicts observation [McNutt, 1983, 1984]. It is therefore worth 
noting that the revised cold plasma input parameters em- 
ployed in our model eliminate this effect here, and consider- 
ing first the current ratio profile for (Sp/M) — 10~'^ profile, 
it is apparent that the hot plasma current is significantly 
larger than that of the centrifugal force over essentially all 
the magnetosphere. However, this is not the case for the 
higher values of (Ep / M), for which the centrifugal force cur- 
rent exceeds the hot plasma pressure current outward of ~54 
and ~34 Rj for {T,*p/M) = IQ-*, and 5 x 10"* mho s kg-\ 
respectively. It should be noted that this does not contrar 
diet the conclusions of McNutt [1983, 1984] and Mauk and 
Krimigis [1987], which were based on Voyager data that 
were obtained at current sheet crossings within 40 Rj and 
that are somewhat sparse beyond ~30 R,j (sec, e.g. Figure 2 
of McNutt [1983]). McNutt [1984] and Mauk and Knmigis 
[1987] computed the ratio of the rotational kinetic energy 
density to magnetic energy density, which can be thought 
of a 'plasma beta for bulk rotation', comparable to the tra- 
ditional plasma beta /3 = (P/Pb), where Pb = B^/2/Uo is 
the magnetic energy density (note they termed this quantity 
M'^, since it is equal to the square of the Alfvcnic Mach mim- 
ber). Achilleos et al. [2010] pointed out that the plasma beta 
for bulk rotation is given by /3cent = (/?cp'^/2/^), and thus 
confirmed that in Caudal's [1986] model the hot plasma beta 
lih dominates the bulk rotation beta /3cent beyond ~40 Rj. 
We have calculated the ratio {pcent / Ph) using our model re- 
sults and, while for clarity we have not plotted the profiles 
in Figure 5, we note that they are very similar to those for 
(jVo cent/jVo h)- AchUlcos ct ttl. [2010] showed that in the 
original model of Caudal [1986], /3cont peaks near ~25 Rj 
at ~16, whereas McNutt [1984] obtained values of ~3 near 
25 R.J. In our results /3cent = 0.75, 6.63, and 24.96 at 25 Rj 
for (Sp/M) = 10^^ 10"'', and 5 x 10"* mho s kg-^ re- 
spectively, with the value for (Ep/M) = 10~* mho s kg~^ 
thus being in most agreement with observations. 

The cold plasma pressure computed in this model is plot- 
ted in Figure 5e, along with power laws fitted by Frank et al. 
[2002] to the pressure values as measured by the Galileo 
spacecraft. These fits are given by 

. _ / l.9 X 10-4 po-^-'^ Pa if Po < 20 Rj 
"^''"^ " \8.6 X 10-« po-i-«^ Pa if Po > 50 Rj , ^ ^ 

shown by the dot-dashed and dashed black lines, respec- 
tively, although we note that in Figure 8 of Frank et al, 
[2002], the scatter in the measured values is generally 
at least an order of magnitude, and in the region 20 < 
(po/Rj) < 50 the points generally lie between the two power 
laws. All three cold plasma pressure profiles are similar 
out to distances of ~15 Rj, beyond which the profiles for 
(Ep/M) = 10~* and 5 x 10~* mho s kg~^ are in best agree- 
ment with the observed profile. 

Figure 5f shows the azimuth-integrated equatorial radial 
current computed from the plasma angular velocity and 
magnetic field profiles using equation 7, where we note that 
for the M-I coupling current equations we explicitly take the 
canonical value of M = 1000 kg s~^, such that Ep = 0.01, 
0.1, and 0.5 mho. The solid lines show the results obtained 
using the magnetic field model discussed here, while the dot- 
dashed lines show the profiles obtained using the empirical 



CAN-KK magnetic field model employed in previous stud- 
ies for comparison. It is evident that the current profiles are 
similar to the results for the CAN-KK model out to ~40- 
50 Rj, beyond which they reduce to smaller values owing 
to the lower values of Fe in the outer region relative to the 
CAN-KK values due to the current sheet outer fringing field. 
Nichols and Cowley [2004] used the midnight sector Galileo 
data of Khurana [2001] to show that the observed values 
of Ip increase rapidly in the inner region, between ~15 and 
25 Rj, before platcauing at ~100 MA at distances beyond, 
out to ~100 Rj (sec, e.g. their Figure 12). It is worth noting 
that the Khurana [2001] data obtained at midnight is un- 
constrained by an assumed magnetopause distance of 85 Rj, 
such that it is not surprising that the decrease in the outer 
region is not evident in those data. This aside, the current 
profile which best fits this pattern is that for Ep = 0.5 mho. 

The resulting field-aligned current at the top of the iono- 
sphere computed using equation 2 is then plotted in Fig- 
ure 5g. In the inner region the currents arc upward and 
peak at similar values to those obtained using the CAN-KK 
field model, at radial distances of ~22, 28, and 33 Rj for 
Ep = 0.01, 0.1, and 0.5 mho, respectively. However, the 
decreasing values of Ip in the outer region result in a re- 
versal of the field-aligned current at ~40-60 Rj, such that 
the current is then downward in the region beyond. Note 
that the oscillation in the Ep = 0.5 mho profile between 
~40-60 Rj is due to the instability in the magnetic field 
model as discussed above. While wc note that such layering 
of upward and downward field-aligned current has been ob- 
served in Jupiter's middle magnetosphere by Mauk and Saur 
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Figure 6. Plot showing (a) the maximum azimuth- 
integrated radial current, (b) the maximum field-aligned 
current density at the top of the ionosphere, (c) the equa- 
torial distance of the maximum field-aligned current, and 
(d) the ionospheric co-latitudes of the maximum field- 
aligned current (joined crosses) and Ganymede footprint 
(joined asterisks), all versus ionospheric Pedersen con- 
ductance. Also shown by the dashed lines are the results 
for a power law current sheet magnetic field. 
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[2007], wo do not wish to infer too much from the structure 
in our results, and simply note that the overall structure 
is that of consistent upward current inward of ~40 Rj and 
downward current outward of ~60 Rj. This confinement 
of the upward field-aligned current to the region inward of 
~40-60 Rj, depending on (Ep/M), is consistent with the 
results of Vogt et al. [2011], who showed using flux equiva- 
lence calculations that the poleward boundary of the main 
auroral oval maps to ~30-60 Rj depending on local time, 
and we also note that Khurana [2001] showed using Galileo 
data that the main oval field-aligned currents fiow inward 
of 30 Rj. Inclusion of local time asymmetry is not possi- 
ble in our axisymmetric model, but the overall results are 
broadly consistent with the observations of Khurana [2001] 
and Vogt et al. [2011]. The downward current in the region 
outward of ~40-60 Rj thus corresponds to the dark polar 
region just poleward of main oval, which typically exists 
on the dawn side but sometimes extends to all local times 
[Grodent et at, 2003b; Nichols et at, 2009]. Note that while 
Nichols and Cowley [2004] showed that the modulation of 
the ionospheric Pedersen conductivity by auroral electron 
precipitation concentrates the peak field-aligned current in 
the ~20-40 Rj region, in their model the field-aligned cur- 
rent was still upward throughout the magnetosphere, albeit 
at low values in the outer region. In addition, while Cow- 
ley et al. [2005] included a region of downward current in 
the outer magnetosphere by design of their specified plasma 
velocity profiles, the results presented here are the first to 
self-consistently produce this downward current region. The 
latter authors also showed that a second sheet of upward 
field-aligned current should exist, associated with the iono- 
spheric flow shear at the boundary between open and closed 
field lines, and indeed it is thought that Saturn's main auro- 
ral oval is due to such a layer between the outer edge of the 
ring current and the open-closed field line boundary [Bad- 
man et al, 2006; Bunce et al., 2008]. Since our model only 
includes closed field lines, we do not consider this second 
layer of upward current and simply note that it will act to 
modify the field-aligned current profiles in the very outer 
region from those computed here. 

Overall, then, it is apparent that the results for 
(Ep/M) = 10"** and 5 x 10"" mho s kg"" provide the best 
agreement with various sets of observations, with the mag- 
netic field (both \Bze\ and B<p) most consistent with the 
latter, and plasma data (angular velocity and pressure) in 
best agreement with the former. It is, however, instructive 
to examine how the M-I coupling currents vary in peak mag- 
nitude and location over a range of values of (Ep /M) , which 
we thus show in Figure 6, again taking here M — 1000 kg s"" 
(note that we consider the effect of changing M in Sec- 
tion 3.2). From top to bottom, the joined crosses in Fig- 
ure 6 show the maximum azinmth-integrated equatorial ra- 
dial current, the maximum upward field-aligned current den- 
sity at the top of the ionosphere, the equatorial radial dis- 
tance of the peak upward field-aligned current, and finally 
the ionospheric co-latitude of the peak upward field-aligned 
current. Also shown for comparison by the dashed lines 
in Figure 6 are results calculated using the analytical solu- 
tion of the Hill-Pontius equation (equation 3) obtained by 
Nichols and Cowley [2003] for a power law current sheet 
magnetic field which maps to a thin latitude band in the 
ionosphere, such as that given by the second term in equa- 
tion 34, thus appropriate for the jovian middle magneto- 
sphere. The analytic power law field result for the maximum 
azimuth-integrated radial current is given by 

Ipmax = 87rEpf2jFo , (38) 

where Fo is the value of Fe at the location of the lat- 
itude band, taken by Nichols and Cowley [2003] to be 
Fo = Fe{70 Rj) ~ 3.22 x 10* nT R^, a representative value 



for the middle magnetosphere current sheet. The maxirrmm 
value for the power law field strictly occurs at pe = oo, 
while at large but finite distances in the rmmcric solution 
using the full empirical field model (e.g. beyond 1000 Rj for 
Ep = 0.5 mho). It is apparent that the mciximum radial 
current computed using the model employed here increases 
less quickly with Ep than for the power law field, i.e. from 
~5 MA at Ep = 0.01 mho to ~242 MA at Ep = 2 mho. 
This occurs since, for a power law current sheet field the 
total azimuthal current increases monotonically toward the 
maximum value given by equation 38 at a rate determined 
solely by the corotation breakdown distance, given by ph 
for a dipole field and by 

for the power law current sheet field. The radial current pro- 
files obtained using the model presented here, however, are 
also constrained by the assumed magnetopause distance and 
thus drop away from the power law profiles at distances in- 
creasingly small relative to ph cs as (Ep/M) increases, such 
that the peak current rises less quickly with Ep than for the 
power law field. 

The maximum field-aligned current current density at the 
top of the ionosphere is plotted in Figure 6b, alongside the 
result for the power law current sheet field, shown by Nichols 
and Cowley [2003] to be 

- 3.05 (^) (^)"" ^hBMj . (40) 

It is evident that the maximum field-aligned current den- 
sity computed here increases with Ep similarly as does 
the result for power law field, i.e. from ~0.004 pA m"'^ at 
Ep = 0.01 mho to ~4.6 pA at Ep = 2 mho, such that 
the latter is a reasonable approximation for the results ob- 
tained here. 

The same is not true, however, for the equatorial radial 
distance of the peak field-aligned current density, shown in 
Figure 6c, in which the distance for the power law field, 

given by 

rises much more quickly than do the results here, which 
increase from 22 Rj at Ep = 0.01 mho to ~44 Rj at 
Ep = 2 mho. This is again due to constraint by the fi- 
nite magnetopause distance in the current sheet field model 
employed here, rather than the power law field which simply 
decreases monotonically toward pe = oo. The ionospheric 
co-latitudes of these peak field-aligned current locations are 
shown in Figure 6d by the joined crosses, along with the 
mapped location of Ganymede's orbit at 15 Rj, shown by 
the joined asterisks. Again, shown by the dashed line for 
comparison is the location of the peak field-aligned current 
for the power law field, given by 
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which indicates that in this case the peak current shifts pole- 
ward as the equatorial radial distance increases, although, as 
shown in Figure 5b, the co- latitude is only weakly dependent 
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Figure 7. As for Figure 5, except with Sp =0.1 mho 
and M = 500 (blue), 1000 (green), and 2000 (red) kg 8"^ 



that in this case the outward transport rate is assumed con- 
stant. Here, we take Ep =0.1 mho, and M = 500, 1000, 
and 2000 kg s~^, typical of the range of values determined 
by various studies [e.g. Hill, 1980; Khurana and Kivelson, 
1993; Delamere and Bagenal, 2003]. Figure 7 thus shows 
the magnetodisc and M-I coupling current system parame- 
ters for constant plasma density in the same format as for 
Figure 5, while Figure 8 shows the results taking the cold 
plasma density to be given by equation 32. It is first evi- 
dent from Figures 7a and 8a that taking Nc oc M acts to 
suppress the divergence of the \Bze\ profiles in the middle 
magnetosphere beyond ~20 Rj. It is, however, just appar- 
ent that for the case with Nc oc M, the higher value of M 
leads to slightly lower equatorial magnetic field strengths 
in the region inside ~40 R.j, i.e. the opposite behaviour to 
the case with constant Nc- This is more evident in Fig- 
ures 7b and 8b, in the former of which the field maps to 
lower co-latitudes for higher mass outflow rates, indicating 
a less stretched held, while in the latter case the held maps 
to higher co-latitudes, indicating a more stretched held. It 
is also worth noting that for the case in 8b the difference 
in held mapping is larger at all radial distances than for 
7b, in which the divergence is only significant outward of 
~15 Rj. This indicates the nature of the centrifugal force 
acting on the plasma in the two cases, which we now discuss. 

Figures 7c and 8c show that the plasma angular veloc- 
ities in the two cases are very similar, with perhaps mod- 
estly increased values in Figure 7c over those in Figure 8c. 
The plasma angular velocity is related to the centrifugal 
force acting on the rotating plasma, which we recall from 
equation 11 is proportional to du^, such that if these two 
parameters are dependent on M, the centrifugal force is pro- 
portional to some power of M, i.e. M'' . An understanding 
of the difference in behaviour between the two cases can 



on the radial distance due to the stretching of the middle 
magnetosphere field lines. However, although the equatorial 
radial distance of the peak field-aligned current in the results 
presented here increases with Ep, above Ep = 0.05 mho, 
the ionospheric co-latitude actually increases slowly with 
Ep, moving from ~16.3° for Ep = 0.01 mho to ~16.8° for 
Ep = 2 mho. This arises since the outward movement of 
peak field-aligned current with increasing Ep is offset by 
the modified mapping of the increasingly stretched magnetic 
field. This can be appreciated by examination of Figures 5b 
and 5g, in which the peak field-aligned current moves out- 
ward for the blue, green and red profiles, respectively, while 
the associated ionospheric mapping profiles also move equa- 
torward, counteracting the outward shift. Considering now 
the co-latitude of the Ganymede footprint, it is evident that 
this is also only very weakly dependent on Ep, moving from 
~17.6° for Ep = 0.01 mho to ~18.0° for Ep = 2 mho. 
This is simply due to the fact that in these runs the mag- 
netic field model is relatively insensitive to changes inside 
~15 Rj, and, as can be seen from Figure 5b, the mapped 
ionospheric co-latitudes of — 15 Rj are very similar for 
all values of Ep. Thus, on the basis of these results, it is 
unlikely that a change of ionospheric conductance is respon- 
sible for the ~3° and ~2° shifts in latitude of the main oval 
and Ganymede footprint, respectively, reported by Grodent 
et al. [2008]. In the following section we therefore examine 
the effect of changing cold plasma number density. 

3.2. Comparison with results taking Nf. proportional 
to M 

We now compare results for which the cold plasma den- 
sity is assumed constant, such that the outward transport 
rate is proportional to M, with results for which the cold 
plasma density is taken to be given by equation 32, such 
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Figure 8. As for Figure 7, except with the cold plasma 
number density given by equation 32. 
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then be obtained if we consider the power law magnetic 
field approximations of Nichols and Cowley [2003]. In this 
approximation, the plasma angular velocity scales with the 
current sheet 'Hill-distance' pn cs given by equation 39, such 
that oj oc M~^^"^. Hence, if the plasma density is indepen- 
dent of the mass outflow rate, we have 7 = — 2/m, while if 
the plasma density is proportional to the mass outflow rate 
we have 7 = 1 — 2/m. Therefore, in the former case 7 < 
for all positive values of m (i.e. for fields which decrease 
in magnitude with distance), such that the centrifugal force 
decreases with increasing M. On the other hand, for the lat- 
ter case we have 7 < for m < 2, such that the centrifugal 
force decreases with increasing M, and < 7 < 1 for m > 2, 
such that the centrifugal force increases with increasing M 
in this case. Thus, examination of Figure 7d, for which 
7 = —2/m, indicates that the centrifugal force is lower for 
increasing M (note that the hot plasma pressure current 
does not differ significantly between the different M cases). 
In Figure 8d, on the other hand, in the inner region where 
the field strength decreases quickly, centrifugal force is larger 
for higher values of M, while in the outer region, where the 
field is very weakly dependent on pe,.the centrifugal force is 
somewhat lower for higher values of M. Physically, the com- 
peting effects of increasing M, i.e. increased plasma density 
but decreased angular velocity, mutually counteract in the 
middle magnetosphere, such that the magnetic field in this 
region becomes relatively insensitive to the value of M. It 
is important to note that the power law approximation is 
not perfectly applicable to the model results obtained hero; 
for example, in the outer region, the field strength increases 
slowly with radial distance, a situation not considered by 
Nichols and Cowley [2003], and for which the power law 
approximations were not designed. Second in the outer re- 
gion, the field does not map to a narrow band in the iono- 
sphere, such that the the approximation conditions do not 
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Figure 9. As for Figure 6, except that here the results 

are plotted versus M, and the points joined by the dot- 
ted lines assume the cold plasma density is independent 
of M, while those joined by the solid lines assume it is 
given by equation 32. 



strictly hold in this region. Hence, while caution should be 
used when comparing with the power law approximation, 
it nevertheless gives a reasonable insight into the behaviour 
of the system. The profiles shown in Figures 7d and 8d 
also indicate why the ionospheric mapping differs between 
the two cases. In the former case, the centrifugal force is 
solely dependent on the plasma angulax velocity, which in- 
side 15 Rj is not particularly sensitive to M, such that in 
this region the azimuthal current profiles, and thus the field 
mapping, do not differ greatly. In the latter case, the cen- 
trifugal force also depends on the plasma density, such that 
the azimuthal current, and thus the field mapping, in the 
inner region varies significantly with M. 

Figures 7e and 8e indicate that taking Nc (x M causes 
the cold plasma pressure to vary more significantly over the 
region inward of ~60 Rj than otherwise. In addition, the 
M = 2000 kg s"^ profile fits the Frank et al. [2002] power 
laws best here, although it should be noted that these pro- 
files are dependent on what reference value for M is used 
and, for example, higher pressure values would be obtained 
if reference values less than 1000 kg s~^ had been taken. 

The azimuth-integrated radial current profiles arc shown 
in Figures 7f and 8f, while the field-aligned current profiles 
are shown in Figures 7g and 8g. Both sets of current pro- 
files in Figures 7 and 8 arc reasonably similar, differing most 
significantly in the degree to which they track the CAN- 
KK results, leading to different peak current values as will 
be discussed further below. As Nichols and Cowley [2003] 
pointed out, for the CAN-KK field model the radial and 
field-aligned currents both tend to values dependent only 
on M in the inner region, and the radial current tends to 
a value dependent on Sp at large distances. In these fig- 
ures the currents thus exhibit 3 distinct profiles in the inner 
region, in contrast to Figures 5f and 5g, although the ra- 
dial current profiles do not converge on a single value in the 
outer region due to the decrease in current intensity owing 
to the increased field strength over the CAN-KK model in 
this region. In both cases the field-aligned current reverses 
from upward to downward between ~45-50 Rj, decreasing 
with radial distance for higher values of M. 

Turning now to Figure 9 we show the magnitudes and 
locations of the peak currents versus M in the same format 
as for Figure 6, except that here the points joined by the 
dotted lines show results taking the cold plasma density to 
be independent of M, while those joined by the solid lines 
show those taking it to be given by equation 32. Figure 9a 
shows that the peak azimuth-integrated radial current in- 
creases with M for both cases, from ~15 MA to ~40 MA, 
and ~5 MA to ^-50 MA for constant Nc and Nc oc M, re- 
spectively as M goes from 100 kg s^^ to 4000 kgs^^, i.e. 
somewhat quicker for Nc oc M. As is evident from Fig- 
ure 7, this arises since, as M increases while Nc is con- 
stant, the magnetic field becomes less stretched due to the 
lower plasma angular velocity, such that the Ip profiles fall 
away from the CAN-KK results at closer distances. Thus, 
the peak currents increase slowly with M. On the other 
hand, for Nc oc M the reverse is true, i.e. the field is more 
stretched for higher M, such that the cm-rent profiles follow 
the CAN-KK results further, increasing the rate at which 
the peak current increases with AI. This behaviour also ac- 
counts for the difference in the field-aligned current profiles 
shown in Figure 9b, in which changes from ~0.2 fiA 
to 0.05 AtA m"^, and --0.01 /jA m"^ to 0.1/iAm"^for 
constant Nc and Nc oc M, respectively as M goes from 
100 kg s"^ to 4000 kg s"^ 

The equatorial distance of the peak field-aligned current 
shown in Figure 9c changes similarly for both cases, i.e. 
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decreasing from ~35 Rj to ~23 R,t, and from ~31 Rj to 
~24 Rj for constant Nc and Nc oc M, respectively as M goes 
from 100 kg s^^ to 4000 kg s^^. However, the difference in 
the field mapping results in the ionospheric co-latitudes plot- 
ted in Figure 9d varying differently with M for the two cases. 
First, as for Figure 6d, for constant Nc the radial motion of 
the peak field-aligned current is offset by the changing field 
mapping such that the ionospheric C9-Iatitude of the peak 
current is only weakly dependent on M, changing from from 
~16.8° for M = 100 kg s"^ to ~16.3° for M = 4000 kg s"\ 
On the other hand, the reverse behaviour of the field map- 
ping for Nc oi M reinforces the radial motion in this case, 
such that the the co-latitude increases more rapidly than for 
the CAN-KK field, moving from ~15.1° for M = 100 kg s~^ 
to ~17.5° for M = 4000 kg s"\ Similarly the co-latitudes of 
the Ganymede footprint change from ~18.0° to ~17.7°, and 
~17.4° to ~18.5° for constant Nc and Nc oc M, respectively 
as M goes from 100 kg s"^ to 4000 kg s~^ Thus, while it 
is difficult to generate the ~3° and ~2° shifts in latitude of 
the main oval and Ganymede footprint reported by Grodent 
et al. [2008], these results suggest that a significant change 
in the iogenic plasma mass outflow rate, combined with an 
associated variation in the cold plasma density in the mag- 
netosphere, possibly as a result of changing volcanic activity 
on lo, is the best candidate for explaining the shift in these 
auroral features. In this case then, the blue image in Fig. 2 
corresponds to an epoch of low volcanic activity, and the red 
image corresponds to an interval of high activity. 

4. Summary and Discussion 

In summary, we have considered the effect of a sclf- 
consistcntly computed magnetosdisc flcld structure on 
the magnetosphere-ionosphere coupling current system at 
Jupiter. Speciflcally, we have incorporated the calculation 
of the plasma angular velocity profile using Hill-Pontius the- 
ory into the model of Caudal [1986], such that the resulting 
magnetosphere-ionosphere currents are computed using val- 
ues of the equatorial magnetic field self-consistent with the 
plasma angular velocity profile. We have thus obtained re- 
sults using a more realistic plasma angular velocity profile 
than that used by Caudal [1986], and we have also included 
updated plasma parameters from Galileo data. We have 
then examined the effect on the system of the values of two 
key magnetosphere-ionosphere coupling current system pa- 
rameters, i.e. the ionospheric Pedersen conductivity Ep and 
iogenic plasma mass outflow rate M . 

We have thus found that the azimuthal current inten- 
sity is dependent on the values of Ep and M . Speciflcally, 
if the cold plasma density is taken to be independent of 
M, we flnd that higher values of the quotient (Ep/M) re- 
sult in a more stretched magnetic field structure with a 
thinner, more intense current sheet due to the increased 
centrifugal force owing to higher plasma angular velocities. 
We thus find that the results for (Ep/M) = ~ 10"" to 
5 X 10~* mho s kg~^ provide the best agreement with vari- 
ous sets of observations, with both \Bzc\ and B.^ data most 
consistent with (Ep/A/) = 5x 10^" mho s kg^^, and plasma 
angular velocity and pressure data in best agreement with 
(Ep/M) = lO^'' mho s kg^^. We discuss a possible reason 
for this discrepancy below. In addition, we find that, while 
the equatorial azimuthal current in the inner region is dom- 
inated by hot plasma pressure, as is generally held to be the 
case at Jupiter, the use of a realistic plasma angular velocity 
profile actually results in the centrifugal current becoming 
dominant in the region beyond ~35-50 Rj, with the exact 
distance depending on the value of (Ep /M) taken. This sit- 
uation similar to that which has been determined for Saturn 
[Achilleos et al., 2010]. 



Overall, the equatorial magnetic field profiles obtained 
are reasonably similar in the inner region to the empirical 
CAN-KK model used in previous studies, such that the cur- 
rents are of the same order as previous solutions obtained 
using this fixed equatorial field strength model. However, 
we show that the outer fringing field of the current disc acts 
to reverse the field-aligned current in the outer region, thus 
reproducing the dark region just poleward of the main oval. 
The confinement of the upward current region to within 
~40-60 Rj is consistent with the recent mapping of Jupiter's 
auroral features to the equatorial plane by Vogt et al. [2011]. 
These authors also determined the location of the open- 
closed field line boundary to be at ~11° co-latitude, a result 
which is also consistent with the ~8-ll° co-latitudes of the 
last closed field line in the model presented here. Further, 
we have found that, while the peak magnitudes of the M-I 
coupling currents arc similar to those which have been de- 
termined previously, we have shown that the location of the 
peak currents differs significantly. For example, the equato- 
rial radial distance of the peak field-aligned current density 
increases with (Ep /M) slower than simply using the CAN- 
KK model. However, if the plasma density is independent of 
M, this outward motion is counteracted by the simultaneous 
stretching of the field, such that the ionospheric co-latitude 
of the peak remains essentially constant. We have there- 
fore also examined the case whereby the plasma density is 
taken to be proportional to M. Hence, we found that in the 
inner region, where the field magnitude decreases quickly 
with distance, the centrifugal force increases with M, while 
in the outer region, where the field in this model does not 
vary greatly with distance, the opposite is true. Overall, 
the competing effects of increasing M, i.e. increased plasma 
density but decreased angular velocity, mutually counteract 
in the middle magnctosphcrc, such that the magnetic field 
in this region is then relatively insensitive to the value of M. 
However, the nature of the centrifugal force in this case is 
such that changes to the field mapping induced by varying 
M now reinforce the associated radial motion of the peak 
field-aligned current, such that the ionospheric co- latitude 
of the peak current varies more significantly, with higher 
values of M corresponding to lower co-latitudes. However, 
very large variations in the plasma mass outflow rate, well 
over an order of magnitude, are still required to reproduce 
shifts comparable to those observed by Grodent et al. [2008]. 

There are various directions in which this work should be 
taken forward. First, the ionospheric Pedersen conductiv- 
ity is assumed for simplicity to be constant, such that feed- 
back effects due to auroral precipitation are neglected. How- 
ever, Nichols and Cowley [2004] showed that precipitation- 
induced enhancements of the Pedersen conductivity signif- 
icantly affect the currents and plasma flows, such that this 
should be taken into account in future developments of this 
model. Similarly, the flcld-aligned voltages required to drive 
the fleld-aligned currents are neglected in the present model. 
The signiflcance of such voltages has been previously de- 
bated [Nichols and Cowley, 2005; Ray et al., 2010], and it 
would be worth determining their effect in the model pre- 
sented here. Third, the Caudal [1986] model neglects the 
effects of plasma pressure anisotropy, which has been shown 
by Paranicas et al. [1991] to be a signiflcant factor in the 
radial stress balance in Jupiter's magnetosphere, and indeed 
has been recently shown to be important at Saturn [Kellett 
et al., 2011]. It is probable that the omission of pressure 
anisotropy in the model is the cause of the discrepancy be- 
tween the plasma and magnetic field data in terms of which 
value of (EJ,/M), i.e. ~ 10"* or 5 x 10"* mho s kg-^ pro- 
vides the best agreement. Development of the model to 
include this effect may produce significant inroads into the 
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problem. A fourth obvious area for further study is to ex- 
amine the effect of the assumed sub-solar magnetopause dis- 
tance. Here we have simply taken the representative value 
of 85 Rj, but observed values range over ~45-100 Rj [Khu- 
rana et al, 2004], and Khurana [2001] presented evidence 
of the solar wind's influence on .Jupiter's magnetic field. 
Caudal [1986] examined the effect of assumed magnetopause 
distance and showed that the larger the assumed distance, 
the more disc-like the field, and thus for different assumed 
magnetopause distances, the field mapping and M-I currents 
will be modified from those presented here, and this should 
be examined in future studies. The location of the mag- 
netopause is governed by the condition of pressure balance 
between a combination of magnetic and plasma pressures on 
one side and shocked solar wind ram pressure on the other, 
with variations typically being taken to be caused by vari- 
ations in the solar wind dynamic pressure (e.g. Huddleston 
et al. [1998]). Cowley et al. [2007] examined the effect of 
solar wind-induced expansions and compressions on the jo- 
vian M-I coupling current system using a prescribed field 
model. In their model a strong compression which reduces 
Rmp from 85 to 45 Rj results in modified field mapping such 
that the peak upward field-aligned current moves poleward 
by ~1°, and the ionospheric mapping of Ganymede's foot- 
print at 15 R.T shifts by ~0.3° Thus, while the solar wind is 
expected to exert some influence on .Jupiter's M-I coupling 
current system, these small latitudinal variations are not 
large enough to account for Grodent et al. 's [2008] observa- 
tions. It seems likely, therefore, that internal factors such 
as the iogcuic plasma disc density are key, but it would be 
illuminating to determine the effect of the solar wind using 
the self-consistent model presented here. Finally, although 
the present model is axisymmetric, the jovian current sheet 
is certainly not [Khurana et al., 2004], such that the effect 
of this should be carefully examined in future studies. 
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